Effect of hydrogel particle additives on water-accessible pore structure of sandy soils: 
A custom pressure plate apparatus and capillary bundle model 
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To probe the effects of hydrogel particle additives on the water-accessible pore structure of sandy 
soils, we introduce a custom pressure plate method in which the volume of water expelled from a wet 
granular packing is measured as a function of applied pressure. Using a capillary bundle model, we 
show that the differential change in retained water per pressure increment is directly related to the 
cumulative cross-sectional area distribution f(r) of the water-accessible pores with radii less than 
r. This is validated by measurements of water expelled from a model sandy soil composed of 2 mm 
diameter glass beads. In particular, the expelled water is found to depend dramatically on sample 
height and that analysis using the capillary bundle model gives the same pore size distribution for 
all samples. The distribution is found to be approximately log-normal, and the total cross-sectional 
area fraction of the accessible pore space is found to be /o = 0.34. We then report on how the 
pore distribution and total water-accessible area fraction are affected by superabsorbent hydrogel 
particle additives, uniformly mixed into a fixed-height sample at varying concentrations. Under 
both fixed volume and free swelling conditions, the total area fraction of water-accessible pore space 
in a packing decreases exponentially as the gel concentration increases. The size distribution of the 
pores is significantly modified by the swollen hydrogel particles, such that large pores are clogged 
while small pores are formed. 

PACS numbers: 



Capillary storage of water is an important property 
that contributes to the plant water availability in soils, 
especially in sandy soils. When the capillary forces are 
strong compared to gravity, the rain water is trapped in- 
side the pores and may be used to support the growth 
of plants. Since capillary forces depend on pore size, the 
amount of capillary water inside a sandy soil is tightly 
linked to its pore structure. An improvement in the wa- 
ter retention of a sandy soil usually couples with a change 
in the water-accessible pore structure. As a popular soil 
additive, superabsorbent hydrogel particles have been 
proven to efficiently enhance water retention of sandy 
soils by swelling and hence locking water inside them- 
selves jTHH] - However, it is yet to be clarified the extent 
to which the improvement is also due to the modifica- 
tion of the water-accessible pore structure caused by the 
presence of the hydrogel particles. 



A standard way to determine the soil water retention is 
to use a pressure plate apparatus, introduced by Richards 
in 1940's [ini E]. The basic idea is to measure the 
amount of solution expelled from a wet soil under a given 
pressure head. For this, a water-saturated soil sample is 
placed on an extraction chamber whose bottom is embed- 
ded with a wet porous plate. When an extra gas pressure 
is applied, the wet porous plate allows soil water to flow 
out but prevents the escape of the compressed gas in 
the chamber. In the past several decades, this apparatus 
has been widely used in soil research to determine the 
so-called soil- water characteristic curve (SWCC), 9{P) 
vs P, which is defined as the ratio of the water volume 
V wa ter retained in soil under a given suction pressure P 



to the initial volume of the dry soil V so u: 

9{P) = K,ater(P) . (1) 

Vsoil 

This is also referred to as the soil-water retention curve 
(SWRC), the degree of saturation, and the volumetric 
water content. The characteristic curve allows a direct 
comparison of the water-holding capacity of soils; see, 
e.g., Ref. |12j for example data and empirical fitting 
forms. The characteristic curve also contains informa- 
tion about the soil pore structure; see, e.g., Ref. [T3] for 
a review of mathematical models. In spite of this body 
of work, the experimental accuracy and reproducibility 
of this approach have been long-standing issues. Studies 
[T3HT5] have shown that the soil water content results ob- 
tained from this technique vary when different operating 
procedures and measuring time scales are used. Most of 
these previous studies focus on high gas pressures, under 
which the real equilibrium state become extremely hard 
to reach and its influence dominates the result accuracy. 
Little attention has been paid to the influence of soil 
packing height, or to operation under low gas pressure. 

Mercury porosimetry [T§H22] is a popular method for 
characterizing the pore structure for rocks, rather than 
soils, but has close parallels to the pressure plate method. 
For mercury porosimetry measurements, the mercury is 
forced to penetrate into the pores of a dry sample as 
the gas pressure is reduced. The variation of the intrud- 
ing mercury volume with the reduced gas pressure gives 
the cumulative pore size distribution of the sample. Data 
analysis relies upon two assumptions: first, that the sam- 
ple pores have a cylindrical geometry; second, that the 
pressure difference due to the sample packing height can 
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be ignored when compared to the applied pressure [2"2T - 
I24j . However, due to the high density of the mercury 
(e.g., 1 cm mercury column corresponds to a pressure of 
1.3 kPa.), the second assumption may cause large devi- 
ations when data is obtained at low pressure values or 
from a relatively high sample packing. Considering the 
similarity between the pressure plate measurement and 
the mercury intrusion measurement, we may convert the 
soil water retention data to the soil pore size distribution 
in the same way, except that the advancing contact angle 
of mercury used in the deduction should be replaced by 
the receding contact angle of water. 

In this paper we build a custom pressure plate appara- 
tus for measuring the volume of the expelled water from 
a soil sample as a function of applied pressure P. Rather 
than use the characteristic curve 0(P) to analyze the re- 
sults, we introduce a new dimensionless parameter, the 
differential expelled water curve 



E(P) 



99 dV w 
A dP ' 



(2) 



where dV w is the incremental volume of water expelled 
by increasing the pressure across the sample from P to 
P + dP, A is the sample cross-sectional area, p is the 
density of the expelled water, and g — 9.8 m/s 2 . To an- 
alyze the results, we develop our own capillary bundle 
model for extracting the area distribution of pore radii 
from the E{P) curves. The bundles are vertical, rather 
than horizontal [53], and their height plays an impor- 
tant role for water retention that, we emphasize, must 
be explicitly accounted for in order to correctly analyze 
pressure plate data. Failure to do so would introduce 
a systematic error, and hence an uncontrolled source of 
irreproducibility, that we make obvious. 

The paper is organized as follows. We begin by describ- 
ing our custom pressure plate apparatus, the model soils, 
and the procedures for taking data. Then we introduce a 
capillary bundle model that directs the extraction of the 
water-accessible pore size distribution from experimen- 
tal data. The validity of the model is demonstrated by 
comparison of results obtained for model sandy soil sam- 
ples with several different packing heights, which have 
extremely different E(P) curves. Finally, the samples of 
uniformly mixed model soils and hydrogel particles are 
examined. The effects of gel concentration, gel size, and 
external confinement on the soil pore structure are de- 
termined respectively. 



I. EXPERIMENT 

Our custom pressure plate apparatus for measuring 
E{P) curves is illustrated schematically in Fig. [I] A 
cylindrical glass column (Knotes, NJ) holds the soil sam- 
ple. It is 30 cm height with a constant inner cross- 
sectional area of A$ = 18 cm 2 , and is designed to safely 
pressurize up to 340 kPa. Two PTFE end fittings with 
20 /im porosity polyethylene bed supports are supplied to 
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FIG. 1: (Color online) Schematic of the custom pressure plate 
measurement set-up. A cylindrical glass column with inner 
cross-sectional area of Ao holds a soil sample packing with 
a height of Ho. It connects to a pump and a burette from 
the bottom. The pump is used to pre-saturate the soil sam- 
ple, while the burette with an inner cross-sectional area of ao 
collects the water expelled out during the pressurization. A 
compressed-gas source is used to pressurize the sample col- 
umn. The amount of the pressure in sample column that is 
higher than the gas phase in burette (AP) is measured by a 
differential pressure sensor combined with a high-resolution 
voltmeter. The pressure across the sample packing is then 
determined as P = AP — pgh w (Eq. |3]l). For the fixed vol- 
ume experiments, centimeter-size balls are added on the top 
of the sample packing to fill the empty space of the sample 
column before pre-saturation. 



seal the top and the bottom of this sample column. The 
bottom of the sample column connects to a gear pump 
(Micropump Inc.), a collection burette(Knotes, NJ), and 
a drain outlet through two three-way valves and Tygon 
tubes. The pump can provide flow rates ranging from 
3 to 60 mL/min. It is used to pre-saturate a soil sam- 
ple by a slow upward infiltration of water from below. 
The collection burette is about 40 cm height with an 
inner cross-sectional area of ao = 7.5 cm 2 . It collects 
the water expelled from the soil sample during the pres- 
surization. The top of the sample column connects to 
atmosphere, a compressed-gas source (cylinder of com- 
pressed N2, Airgas Inc.), and a differential pressure sen- 
sor (26PCA, Omega, CT) through two three-way valves. 
The output of the pressure sensor (AP) is measured by 
a voltmeter (Keithley Inc.) with a resolution of 0.1 mV. 
We calibrate the pressure sensor by water columns with 
controlled heights and get a linear dependence with a 
sensitivity of 0.29 kPa/mV. 
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A. Materials 

The model sandy soil we choose are glass beads (Pot- 
ters Industries, PA) with a diameter of 2 mm (±10%). 
To clean the surface, they were first burnt in a furnace at 
500° C for 72 hours and then soaked in a 1 M HC1 bath for 
an hour. After that, the beads were rinsed with DI water 
thoroughly and baked in a vacuum oven at about 110°C 
for 24 hours. The dry glass beads have very hydrophilic 
surfaces. 

The hydrogel particles used in the experiments 
are a commercial product provided by Degussa 
Inc.(Stockosorb SW), made by grinding a bulk gel. The 
particle shape is randomly faceted but compact. The 
main chemical component of these particles is cross- 
linked polyacrylamide-co-potassium acrylate. A small 
amount of salts is present from the industrial polymer- 
ization. If allowed to freely swell in DI water (the salt 
concentration in the final fluid is less than 10~ 3 M) under 
atmosphere, a gel particle can hold several hundred times 
of water than its weight in dry. And 95% of the absorbed 
water is available to plants |25) . In our experiments, two 
different sizes of dry hydrogel particles are chosen. They 
come from the same sample bag but were sieved by differ- 
ent sized copper meshes. The smaller ones (0.2 — 0.3 mm 
diameter) are used in most of the measurements, while 
the larger ones (0.9 — 1.1 mm diameter) are only used 
for comparison. For comparison, the tetrahedral hole for 
our glass beads has diameter 0.225 x 2 mm = 0.45 mm. 

Since water is pulled upwards against gravity into a 
dry hydrophilic sample, other important parameters for 
our system include the liquid mass density p — 1 g/cm 3 , 
the liquid-gas surface tension 7 = 73 dyne/cm, and the 
contact angle 8 = 0° between the liquid-gas interface and 
the hydrophilic grain surfaces. 



B. Procedures and example data 

In this subsection we describe how our samples are 
prepared and how our apparatus is used to measure both 
the pressure P and the dimensionless differential expelled 
water curve E(P) vs P. 

To prepare a soil packing, we first clean and dry the 
sample column to make sure that the inner surface is 
hydrophilic. A piece of filter paper (Whatman, NJ) is 
added on the bottom and the dry granular sample is 
poured into the sample column carefully, 1-2 cm height 
each time, until reaching the desired packing height Hq. 
The sample column is slightly patted every time to form a 
random close packing. Two different conditions are used 
for each mixture: fixed volume and free swelling. For 
fixed volume, centimeter-size plastic balls are added on 
the top of the granular packing to fill the empty space in 
the sample column and maintain the packing at constant 
volume against the pressure of the swelling hydrogel par- 
ticles. For free swelling, no plastic balls are added and 



the sample column has enough empty space for the mix- 
ture packing to freely expand. 

The packing is then pre-saturated with DI water from 
the bottom at a slow flow rate of 3 mL/min with the 
top open to atmosphere. For the pure glass beads pack- 
ings, the pre-saturation procedure does not modify the 
pore structure and the extra confinement has no effect on 
the results. However, for the mixed packings of hydro- 
gel particles and glass beads, the swelling hydrogel parti- 
cles tend to expand the packing during the pre-saturation 
procedure, so free versus fixed volume conditions are dif- 
ferent. After pre-saturation, each mixture packing is left 
for 24 hours to ensure the full swelling of the hydrogel 
particles. 

The final preparation step is to open the top of the 
sample and collection columns to atmosphere, and to 
drain liquid in the collection column down to the same 
level as membrane filter at the bottom of sample column 
- i.e. the location indicated by a horizontal line labeled 
Po in Fig. [l] Once this level is reached, the drain outlet 
is closed; however, over the course of several hours, the 
liquid level in the sample column falls and the level in the 
collection column rises. This extra liquid is then drained, 
and the whole process is iterated as many times as nec- 
essary so that the liquid level in the collection column 
remains constant and even with the bottom of the sam- 
ple column. When this is finally achieved, both columns 
are sealed at top. Referring to the quantities labeled in 
Fig. [l] the initial conditions are thus such that Po equals 
atmospheric pressure, P = 0, AP = 0, h w = 0, and h 
equals the equilibrium capillary rise of water pulled into 
the sample against gravity. 

The dimensionless expelled water curve, E(P) vs P, is 
now measured by pressurizing the space above the sam- 
ple and measuring the resulting expulsion of water from 
the increase in the height h w of water in the collection 
burette. This is done in a step-wise fashion, by repeat- 
edly bleeding in a small quantity of compressed gas (N 2 ) 
and then waiting for the liquid levels to come to equilib- 
rium before h w is recorded. The two quantities directly 
measured are thus the pressured difference AP between 
the columns, and the height h w of liquid in the collec- 
tion burette. In order to deduce the pressure P across 
the sample, first note that the pressure in the collection 
burette at the level of the membrane filter is the same 
value, which we call Po, as at the bottom of the sample. 
This reference pressure is now greater than atmospheric, 
but its value is not of interest. The gas pressure in sam- 
ple column equals Po -I- P, and the gas pressure in the 
collection burette equals Po — pgh w (see Fig. [I]). The gas 
pressure in the sample column is also greater than that 
in the collection burette by the measured quantity AP. 
Altogether this gives the pressure across the sample as 

P = AP-pgh w . (3) 

Note that P gives the amount by which the pressure is 
greater in the gas above, than in the liquid underneath, 
the sample - and hence is sometimes referred to as "sue- 
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tion". In order to deduce E(P), note that the incremen- 
tal volume dV m of water expelled by a small increase in P 
is simply the product a dh w of the inner area of the col- 
lection burette times the change in collected water level. 
The dimensionless differential expelled water curve for 
our apparatus is thus 



E{P) = (a /A )pg dh w /dP, 



(4) 



which may be found by numerical differentiation of h w 
versus P data. 

Example data for h w versus P obtained by the above 
procedures are shown in Fig. [2^,, for 2 mm glass beads 
(no hydrogel particles) packed to different heights Hq as 
indicated in the legend. For all, h w increases monotoni- 
cally with P towards an asymptotic value corresponding 
to the complete expulsion of all water from the sample. 
To reach this limit, the presence of a wet filter paper 
membrane beneath the sample was necessary to prevent 
the penetration of compressed gas out from underneath 
the sample. Note that the results for the highest packing 
heights are nearly identical, and display an initial rise 
from zero that is linear. But the results for the lowest 
packing heights are sigmoidal in shape, and asymptote 
to values that decrease for smaller H . Thus the soil- 
water characteristic curves (SWCC), given by Eq. as 
0(P) — I — h w (P)/h w (oo), clearly depend on packing 
height. This can be understood as follows. When the 
packing is taller than the equilibrium capillary rise, the 
upper portion of the sample is dry and has no influence 
on water retention. When the packing is smaller than 
the equilibrium capillary rise, the entire sample is wet 
and a minimum height-dependent pressure head must be 
exceeded in order for water expulsion to commence. This 
intuition will be made quantitative with a capillary bun- 
dle model in the next section. 

The differential dimensionless expelled water curves, 
E(P) versus P, obtained from the example h w versus P 
data and Eq. Q, are plotted in Fig. ^p. For numerical 
differentiation, we first smooth the data using a LOESS 
algorithm in Igor with a smoothing parameter carefully 
chosen so that the smoothed curves do not systemati- 
cally deviate from the raw data. Then we use central 
finite differencing. For small packing heights, the E(P) 
curves rise to a peak and then fall toward zero as P in- 
creases. For the two tallest packings, the E(P) results 
are indistinguishable, monotonically decreasing with P, 
are equal to a constant fo ~ 0.34 for small P. The value 
of fg represents the total cross-sectional area fraction of 
the water accessible pore space, as will be shown next 
with a capillary bundle model. 



II. CAPILLARY BUNDLE MODEL 

To extract physical meaning from the differential E(P) 
versus P water expulsion curves, we now construct a 
model in which the granular packing is pictured as a set 
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FIG. 2: (Color online) (a) Height of water in the collection 
burette, and (b) dimensionless differential water expulsion 
parameter, versus applied pressure for 2 mm diameter glass 
beads packed to different heights Ho as labeled. As shown in 
Fig. [T] Ao and ao are the inner cross-sectional areas of the 
sample column and the burette, respectively The y-axis in 
(b) represents the expelled water per pressure increment per 
sample area, made dimensionless by the fluid density p and 
g = 9.8 m/s 2 . In (a) the solid curves represent the smoothed 
data used for differentiation. In (b) the solid curves are a 



simultaneous fit to Eq. (Ill for a log-normal pore size distri- 



bution with f(r) given by Eq. (15 1, and the dashed curves 
are the expectation for the low packing height data based on 
the high packing height data and Eq. (111. The plateau in 



(b) labeled fo represents the cross-sectional area fraction of 
water-accessible pore space in the packing. 



of vertical capillary tubes of height Hq and with some dis- 
tribution of radii r, as depicted schematically in Fig. [3] 
Such "capillary bundle" approximations may seem rather 
severe and uncontrolled, but they have a long history of 
use in the modeling of fluids in porous media HH 
[29] . Here, the rise h of liquid into a tube of radius r may 
be computed by considering how the pressure increases 
from Pq in the liquid at the bottom of the sample to 
Pq + P in the gas above. In going upwards from the bot- 
tom to a height h just below the liquid-gas interface, the 
pressure drops according to Pascal's law by pgh. And in 
crossing the interface, the pressure goes up according to 
Laplace's law by where 7 is the liquid-gas surface 
tension and where complete wetting is assumed. In other 
words, the gas pressure Pq + P above the sample is equal 
to Pq — pgh + 27/r, and this gives the capillary rise as 



2j_ 
pgr 



P 

pg' 



(5) 



The first term represents the usual capillary rise formula, 
which would be multiplied by the cosine of the contact 
angle for the case of partial wetting. The second term 
represents the reduction in height due to an applied pres- 
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FIG. 3: (Color online) Schematic of the capillary bundle 
model. In the model, a bundle of capillaries with height of Hq 
is placed vertically in a water reservoir. The pressure in the 
water reservoir is Po and the pressure in gas phase is Po + P 
(P > 0). Here, r sma ii = 2j/(P + pgHo) sets the radius limit 
for the filled capillaries; ri arge — 27/P sets the radius limit 
for the empty capillaries; and the capillaries with radius be- 
tween them are the active ones. Only the water in the active 
capillaries is expelled out when an small pressure increment 
dP is applied to the system. 



are active, the cross-sectional area fraction is 
fo = {ao/A )pg dh w /dP, 



(8) 



which is recognized as our dimensionless differential wa- 
ter expulsion parameter. This holds for tall samples and 
low pressures, for which the E(P) versus P sample data 
in Fig. are indeed constant. For that sample, the in- 
ferred total cross-sectional area fraction of the pore space 
may be read off the graph as fo = 0.34 ± 0.01. Note that 
this argument relies only upon continuity and the second 
term in Eq. ([s]); therefore, we believe its validity tran- 
scends any limitations of the capillary bundle approxi- 
mation. 

Now we generalize to the case that only some of the 
tubes are active. For this we introduce a new quantity, 
the fraction f(r) of the cross-sectional area having pores 
with radii less than r. By definition, /(r) increases mono- 
tonically from and asymptotes to fo as r increases from 
zero to infinity. Also by definition, /(r) is a cumulative 
distribution function and therefore the associated prob- 
ability distribution function is 



PDF(r) 



dr ' 



(9) 



sure (or suction), and is independent of r and wetting 
properties. 

Note that Eq. ^ holds only if r is neither too small 
nor too large. In particular if r is smaller than 



' small 



2" 



P + pgH 



(6) 



then the tube is too short and the rise of fluid will be 
pegged at h = H , with water entirely filling the tube. 
And if r is larger than 



n a 



27 
P' 



(7) 



then the applied pressure is too great and the rise will be 
pegged at h = 0, with all water completely expelled from 
the tube. Tubes of radii in the range r sma u < r < ri arge 
are partially filled with fluid, and are "active" in the sense 
that their filling height h responds to pressure changes 
according to Eq. fj5]). 

The simplest case to consider, first, is a sample with 
tall packing height and with small applied pressure - 
then all tubes are active. Therefore, in response to a 
small increase dP of applied pressure, the change in rise 
height for all tubes is given by differentiating Eq. (J5j) as 
dh = ~dP/(pg). The resulting volume of expelled wa- 
ter is — foA dh where A is the cross-sectional area of 
the sample and fo is the total cross-sectional area frac- 
tion of the pore space. By continuity, the change dh w in 
height of liquid in the collection burette is such that the 
expelled volume equals aodh w . In other words, we have 
aodh w = — foAo[—dP/(pg)]; therefore, when all tubes 



which is normalized to fo rather than to one. As above, 
the volume of water expelled by an increase dP of applied 
pressure is aodh w and equals the active area times the de- 
crease in liquid level inside the sample, — dh = dP/(pg). 
Whereas before the active area was foAo, it is now more 
generally [f{ri arge ) - f(r sma ii)]A where the term in 
square brackets is the fraction of active area with radii 
in the range r smaU < r < n arge that obey Eq. Alto- 
gether the capillary bundle model thus gives 



(«oMo)P5 



dh w 
dP 



f (riarge) ~ f(r s ), (10) 
27 



J ' P J / VP+ W //n 



, (11) 



2-f/(P+pgH ) 



2 7 /P 



PDF(r)dr. (12) 



where the left hand side is recognized as our dimension- 
less differential water expulsion parameter, E(P), and 
the large and small radii are taken from Eqs. (|6][7]). Note 
that the right hand sides reduce to fo, and Eq. ([8]) is 
recovered, in the limit that all pores are active. These 
equivalent expressions are the main result of the capil- 
lary bundle model. In essence, the raw data from pres- 
sure plate measurements of h w versus P are seen to be 
directly linked to a double integral of the cross-sectional 
area distribution of the pore radii. 

As a remark, note that while our model thus shows 
that the cross-sectional area distribution of pore radii is 
the key structural quantity accessible from water reten- 
tion/expulsion data, prior work has been in terms of the 
volumetric distribution of pore radii [13J. 
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III. MODEL SANDY SOIL 

In this section we describe how to use the capillary 
bundle model, Eqs. ( To][l2 ), to deduce pore size infor- 
mation from experimental data for E(P) versus P. And 
we demonstrate the procedures using the example data 
of Fig. [5] for packings of 2 mm glass spheres of various 
heights. 



A. Direct fitting 

One straightforward method of analysis is to assume 
a particular form for /(r) and then simply fit E{P) 



(an / 'Aq) pgdh w I dP data to the right hand side of Eq. (111. 



Here we try three different forms, for which the PDF is 
a peaked function: 



f(r) 

m 

fir) 



/ [1 - e-W**> a ], 
/o[l-(l + (r/r f)e-^\ 
fa 



V2 



x 1 e 



da;. 



(13) 
(14) 

(15) 



All three cumulative distribution functions are sigmoidal 
in shape, rise from /(0) = 0, and asymptote to /n as 
r — > 00. The last form, Eq. (15), is the cumulative dis- 
tribution function for a log-normal PDF. These forms 
are used to directly, and simultaneously, fit the data for 
all five packing heights in Fig. [2}d by keeping f = 0.34 
fixed and adjusting the other two parameters. The fits 
are all satisfactory and give ro = 0.40 mm and a = 4.2 
for Eq. ([13]); r = 0.30 mm and j3 = 2.7 for Eq. ([llf; 
and ro = 0.36 mm and Co = 0.29 for Eq. (15). The log- 
normal fits are displayed in Fig. [2]d and have the small- 
est chi-squared deviation of the three candidate forms. 
A log-normal form is consistent with simulation results 
for the pore size distribution for a random packing of 
spheres [30l432j . We emphasize that, while the E(P) ver- 
sus P curves are all very different and height dependent, 
in effect they all give the same pore radii distribution. 
The good agreement of the simultaneous fits for all data 
sets demonstrated both the consistency of our data and 
validity of the capillary bundle model. 



B. Extraction of f(r) for tall packings 

For large enough packing heights Hq, such that the 
capillary rise of liquid never extends to the top of the 
sample, f(r sma u) vanishes and the right hand side of 
Eq. (11) reduces to f(2j/P). Then the cumulative area 



fraction of pore space with radii less than r is given by 



fir) 



2-f/P 



(ao/Ao)pg dh w /dP, 



(16) 



and a plot of f(r) versus r is obtained directly by plotting 
E(P) = (ao/A )pgdh w /dP data versus 2j/P. Results 
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FIG. 4: (Color online) Variation of (a) cumulative area frac- 
tion of water-accessible pore space with pore radii smaller 
than r and (b) associated probability density, which is nor- 
malized to the total area fraction /o of the pore space. These 
results are based on the water expulsion data of Fig. [2] for 
2 mm glass beads with packing heights Ho as labeled. These 
samples are tal l en ough that f(r) could be extracted directly 
by use of Eq. (16 1. The solid curve is the log- normal form, 
Eq. ( |15[), found by simultaneous fit of Eq. (Ill to all data 
in Fi g. |2| T he dashed curves are similarly obtained fits to 
Eqs. ( 13|14[). The PDF equals d//dr, and the data points 



were computed by finite differencing with no further smooth- 
ing. 



are shown in Fig. [3Ji based on the Fig. [2]d water expul- 
sion curves for the two tallest packings. For comparison, 
the direct fitting results for the three candidate sigmoidal 
forms are also included. The associated probability dis- 
tributions obtained by differentiation, for both the data 
and the fits, are shown in Fig. Hh. 



C. Extraction of PDF(r) for short packings 

For small enough packing heights, such that pgHo <C 
P, the right hand side of Eq. ( 10 ) may be well approx- 
imated as d//dr = PDF(r) evaluated at r ~ 2j/P and 
multiplied by r large - r smaH w 2jpgH /P 2 . Thus the 
radius distribution for the area fraction of pore space is 
given as 



PDF(r) 



2-y/P 



(a /A )pg dh w /dP 
(2 1 /P)(pgH Q /P) • 



(17) 



In other words, a plot of PDF(r) versus r is ob- 
tained directly by dividing the E(P) data by the length 
(2-//P)(pgH /P) and plotting versus 2-y/P. The right 
hand side of Eq. (17) is computed for the three shortest 
packing height data of Fig.[2j for the 2 mm glass spheres, 
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FIG. 5: (Color online) The distribution of the area fraction of 
pore space for 2 mm glass beads, deduced from Eq. ( 17 1 using 
the water expulsion data of Fig. [2] at the three shortest pack- 
ing heights Ho as labeled. The log-normal fit was obtained 
previously, by the simultaneous fit of all data in Fig. [5] [Note: 
the data cut-off used in this figure is set to be P = 0.8 kPa] 



both the concentration and the size of additives. Small 
gel particles (0.2 — 0.3 mm axis in dry) are used to ex- 
amine the influence of gel concentration. In dry, they 
can fit into the existing soil pores without disturbing the 
soil matrix, since the tetrahedral hole for 2 mm diam- 
eter beads is about 0.45 mm. After pre-saturation of 
the sample, however, their maximum swelling size may 
exceed the pore size. Whether they can reach this size 
or not depends on the strength of the soil matrix con- 
finement during pre-saturation. Fig. [6] shows the E(P) 
data at four different gel concentrations in both (a) free 
swelling and (b) fixed volume conditions. The data for a 
"no gel" packing is included for comparison. As the gel 
concentration increases, it become harder and harder to 
expel water out of the soil packing. Note that there is no 
dramatic difference between free swelling data and fixed 
volume data. 



and is plotted in Fig. [5] along with the log-normal dis- 
tribution found from previous fits. We see that results 
become spurious at small r. But more importantly, for 
larger r, we see that the results underestimate the ex- 
pectation but become progressively better for the smaller 
packing heights. Therefore, analysis of water expulsion 
curves with Eq. (17) would requires even shorter samples 
than measured here. 



IV. SOFT HYDROGEL PARTICLES AS SOIL 
ADDITIVES 



A. Gel concentration 

When hydrogel particles are uniformly mixed into 
sandy soils, the pore structure is modified according to 



Direct fitting is applied to extract the water-accessible 
pore structure in these mixture packings. Data in 
Fig. |6^a) and in Fig. [6j^b) are fitted simultaneously by 
following function respectively: 



f(r) = fo 
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(18) 



This is a combination of the cumulative distribution func- 
tions for two log-normal PDFs. The first one represents 
the pores existing between soil particles; while the second 
one represents the pores that may exist between gel par- 
ticles or between gel particle and soil particle. In each fit, 
we fixed ro = 0.36 mm and <7o = 0.29 from the previous 
section on pure glass bead packings, but let the values of 
r\ and <j\ vary simultaneously for all four set of mixture 
packing data. The parameter S represents the percentage 
of the unoccupied soil pores in a mixture packing thus its 
value depends on gel concentration. It is allowed to ad- 
just freely for each set of data, fo is the maximum value 



of f(r) at r — > oo and represents the total area fraction 
of water-accessible pore space in a packing. It also is 
allowed to adjust freely in fits to each set of data. 

The fitting curves are shown in Fig. [6] by solid curves. 
We determine that the value of r\ is 0.18 mm in free 
swelling condition and 0.19 mm in fixed volume condi- 
tion. The value of <j\ is 0.42 and 0.39 for free swelling 
condition and fixed volume condition respectively. The 
results obtained in different conditions are very close to 
each other. The pores existing between gel particles or 
between gel particle and soil particle are only about half 
of the size of the major pores in soil matrix. For fur- 




0.0 0.5 1.0 1.5 2.0 

P (kPa) 



FIG. 6: (Color online) Dimensionless differential water expul- 
sion parameter, versus applied pressure for 2 mm diameter 
glass beads with different concentration small hydrogel par- 
ticles in (a) free swelling condition and in (b) fixed volume 
condition. The solid curves are a simultaneous fit to Eq. ( | 1 1 [ ) 
for a combined log-normal pore size distribution with f(r) 
given by Eq. (18) . The fitting parameters are shown in Fig. [7] 
and Fig. [8] 

ther discussion, we plot the PDF(r) functions for packing 
with different gel concentrations at different conditions in 
Fig. |7ja) and (b) respectively. We clearly see that as the 
gel concentration increases the height of the major peak 
decreases monotonously and finally disappears when gel 
concentration exceeds 0.1 wt%. This corresponds to a gel 
to bead number ratio of about 2:3. We also notice that 
the height change of the secondary peak does not follow 
the same trend as the major one. 

The behavior as a function of gel concentration is sum- 
marized in Fig. [8j Fig. [8]ja) shows the packing height, 
Ho, the results of which are used in the fits. In dry, 
all the mixture packings have the same height. When 
pre-saturated in free swelling condition, the height of the 
mixture packing varies as the gel concentration varies 
by up to about thirty percent. When pre-saturated in 
fixed volume conditions, all concentration packings are 
of course forced to maintain their original height. 

Fig. [8^b) shows the variation of the total area fraction 
of pore space with gel concentration. In both conditions, 
the total area fraction /q decays exponentially as gel con- 
centration increases. This result is consistent with our 
prior studies on the water permeability in the same mix- 
ture system |33j . Again, no significant difference is seen 
between two conditions when correct packing heights are 
applied. The characteristic value of gel concentration in 
the exponential fits is around 0.1wt%, which corresponds 
to the loss of the primary peak in the PDF seen in Fig. [7] 

Fig. [8^c) shows the variation of S occupy , the effective 
percentage of soil pores blocked by gel particles. This is 




r (mm) 



FIG. 7: (Color online) The distribution of the area fraction of 
pore space for 2 mm glass beads with different concentration 
small hydrogel particles in (a) free swelling condition and in 
(b) fixed volume condition, deduced from the simultaneous 
fit of all data in Fig. [6] In both conditions, the major peak, 
whose location and width are fixed to be the result in the 
model sandy soil alone (ro = 0.36 mm and <to = 0.29), reduces 
its height as gel concentration increases; while a secondary 
peak grows in small r region at the same time. In free swelling 
condition, the secondary peak is located at r 1 = 0.18 mm with 
o i = 0.42; in fixed volume condition, it is at n = 0.19 mm 
with ai = 0.39. 

defined as 

S occupy = 100 (l - , (19) 

where 0.34 is the value of /o at c ge i — 0. From the figure, 
we see that the value of S occupy grows linearly, until all 
soil pores are filled, and thus can be fit to an exponential 
function of bead concentration. The gel concentration 
required to fill all soil pores is around 0.1wt%, consistent 
with the characteristic value obtained from the exponen- 
tial fit in Fig. [8jb) . When the gel concentration exceeds 
this value in free swelling conditions, the swollen gel par- 
ticles may occupy extra space by expanding the packing. 
But in fixed volume conditions, the swollen gel particles 
have to share that soil pores with others and cannot swell 
to their desired size. 



B. Gel size 

Lastly, to probe the influence of gel particle size, we 
fixed the gel concentration to be 0.1 wt% and compare 
the results for large gel particle additives (0.9 — 1.1 mm 
axis in dry) to those above for the small gel particle ad- 
ditives (0.2 — 0.3 mm axis in dry) in the same external 
conditions. The large gel particles have a dry size larger 



9 



10 



E 

o 



„° 4 




0.4 



0.3 

0.2 

0.1 

0.0 
100 

80 

60 

40 

20 



J 




■ 


7 


□ 

□ 

□ 

o o o 




□ 


Free swelling 




o 


Fixed volume 


(a) ~: 


1 







1 


1 1 1 I I 1 1 




: ^ 


5%^, /„ =0.02 + 0.32e" 1 "' /c " 3 


': 


1 


= 0.03 + 0.31<"~ r ""'°" 


(b) ; 


i 

r ^ 


= 100(1 - e - c ' l,m ^) == , p 


-. 




Smm =100(1- <.-'■*""*») 








(c) 1 









0.00 



0.05 



0.10 



0.15 



0.20 



0.25 



FIG. 8: (Color online) (a) Packing height and parameters de- 
termined from the simultaneous fit shown in Fig. [6] (b) the 
total area fraction, and (c) effective soil pore occupation rate 
by gel particles. The error bars are smaller than the symbol 
sizes. As gel concentration increases, the total area fraction of 
the pore space decays exponentially for both free swelling con- 
dition and fixed volume condition. The soil pore occupation 



rate is defined by Eq. (19 1. It grows almost linearly as gel con- 
centration increases until all major pores are fully-occupied, 
and hence can be fit to an exponential. 



than the average soil pore size, and thus may perturb the 
soil matrix structure even in dry. After pre-saturation, 
the free swelling packing with large gel particles expands 
to a height of Ho = 7.4 cm while the fixed volume packing 
maintains a height of Ho = 5.1 cm. 

Fig. |9^a) shows the E(P) data for packings with dif- 
ferent size gel particles in free swelling and fixed volume 
conditions, along with the "no gel" data for comparison. 
Under fixed volume, note that the curves are very similar 
for the two different gel particle sizes. Thus, the parti- 
cle size is relatively unimportant, presumably as long as 
it is not very much greater than the bead size. Instead, 
the weight percent of the additives is the more important 
parameter for affecting behavior. 

Fits to the bimodal log-normal pore size distribution 
of Eq. (18 1 are also included in Fig. [9]ja). Since there is 
only one set of data in each condition, we fix the values of 
rn, er , Tx, and <j\ to be the same as those obtained from 
the fits of small gel packings. In free swelling condition, 
the total area fraction is determined as /n = 0.22 and the 

= 69%. 



occupy 



soil pore occupation rate is determined as S, 
Comparing to the corresponding small gel packing, the 
large gel one has a higher total pore area fraction of space 
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FIG. 9: (Color online) (a) Dimensionless differential water 
expulsion parameter versus applied pressure, and (b) the dis- 
tribution of the area fraction of pore space for 2 mm glass 
beads with 0.1 wt% large and small hydrogel particles in both 
free swelling condition and fixed volume condition. The solid 
curves are a fit to Eq. |TT| for a combined log-normal pore 
size distribution with /(r) given by Eq. (18 1. The fitting pa- 
rameters for packings with small gel particle additives have 
been shown in Fig. [5] For packings with large gel particle 
additives, the values of ro, co> n, and o\ are fixed to be the 
same as those obtained from corresponding small gel pack- 
ing. In free swelling condition, the packing height is mea- 
sured as Hq = 7.4 cm, the total area fraction is determined 
as fo = 0.22, and the soil pore occupation rate is determined 



as Soc 



= 69 %. In fixed volume condition, the packing 



height is measured as Ho = 5.1 cm, the total area fraction 
is determined fo = 0.15, and the soil pore occupation rate is 

S occupy = 74%. 



and a lower soil pore occupation rate. In fixed volume 
condition, we determine that the total area fraction is 
fo = 0.15, which is very close to that of a small gel 
packing. The soil pore occupation rate is determined 
as S 'occupy = 74%, which is also lower than that in the 
small gel packing. 

The resulting fitting function for the PDFs are plotted 
in Fig. |9jb). From this figure, we clearly see that under 
free swelling condition more soil pores are left unoccupied 
in the packing when its small gel particle additives are 
replaced by the same mass of large gel particles. How- 
ever, this is a relatively small effect. By contrast, the 
role of external confinement is more important. Namely, 
a strong external confinement helps the gel particles to 
swell into and efficiently block the pore spaces. 
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V. CONCLUSION 

In summary, we developed a custom pressure plate ap- 
paratus for measuring the expulsion of water from model 
soils as a function of applied pressure, and we developed a 
capillary bundle model to analyze the results in terms of 
pore-scale structure. We verified the apparatus and the 
model by obtaining consistent results for the pore struc- 
ture for packings of different heights, where the soil- water 
characteristic curves are all different. And we applied our 
methods to study the effect of superabsorbent gel particle 
additives, which can swell to block the pores. 

One general conclusion is that the height of the sample 
can strongly affect experimental results. Only for very 
tall samples, not initially filled completely via capillarity, 
are data independent of sample height. This point is not 
widely appreciated in prior experiments, where sample 
heights tend to be small and are often not even reported, 
much less varied. This point is also not widely appre- 
ciated in prior modeling efforts, none of which to our 
knowledge accounts for the influence of sample height. 

The capillary bundle model we developed here has gen- 
eral significance for several reasons. First, the sample 
height is now explicitly included in the analysis. Sec- 
ond, it points to the importance of two key concepts 
not considered in prior work: (1) The differential dimen- 
sionless expelled water curve, E{P) versus P, defined 
by Eq. |2]) as the incremental volume dV w of water ex- 
pelled per pressure increment dP, made dimensionless 
by the sample cross sectional area Aq, the density p of 
water, and g — 9.8 m/s 2 ; and (2) The cumulative cross- 
sectional area fraction f(r) of pores with radii less than 
r. Whereas prior experiments focus on the volumetric 
water content and, prior modeling is in terms of the vol- 
umetric distribution of pore radii, the capillary bundle 
model shows how E(P) is simply and directly related to 
f(r) via 

where 7 is the liquid-air surface tension and Hq is the 



sample height. In words, the dimensionless differential 
volume of water expulsion per pressure increment is equal 
to the cross-sectional area fraction of active pores. This 
expression makes explicit both the connection between 
pressure plate data and pore-scale structure, as well as 
the role of sample height. 

In terms of our specific model-soil systems, we can 
make several conclusions based on our new data and 
its analysis via the capillary bundle model. For ran- 
domly packed monodisperse spherical beads, the area 
distribution of the pore radii is well-approximated as a 
log-normal, peaking around one third the sphere radius 
and extending with full-width half-max between about 
0.2 — 0.5 sphere radii; the total cross-sectional area frac- 
tion of the pore space is found to be 0.34. With the ad- 
dition of superabsorbent hydrogel particle additives, this 
pore space is modified by the swelling of the gels. The 
particle size is not crucial, provided it is not very much 
larger than the bead size. Then only about 0.1 weight 
percent is required to clog up the original pores between 
beads and reduce the water-accessible pore space. This 
can have a very significant effect on the permeability of 
the medium |33j . 

We hope that the general experimental and theoretical 
methods presented here may find future use for studying 
both natural soils, as well as the model systems consid- 
ered here. And we hope our work provides a general set 
of tools for reliable quantitative testing of trial additives 
designed to improve water usage efficiency in agriculture. 
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